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Abstract 

Yang's r] pairing operator is generalized to explore off-diagonal long-range 
order in the Hubbard bilayer with an arbitrary chemical potential. With this 
operator and a constraint condition on annihilation and creation operators, 
we construct explicitly eigenstates which possess simultaneously three kinds 
of off-diagonal long-range order, i.e., the intralayer one and the interlayer one 
for on-site pairing, and that for interlayer nearest-neighbor pairing. As in 
the simple Hubbard model there is also an SO (4) symmetry, with the the 
generators properly defined. A sufficient condition leads to at least one of 
the above three kinds of off-diagonal long-range order. A constraint relation 
among different kinds of off-diagonal long-range order is also given. There 
exists a triplet of collective modes if the U(l) symmetry of a subgroup is 
spontaneously broken. 
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I. Introduction 



Off-diagonal long-range order (ODLRO) is essential for the phenomena of super- 
conductivity and superfluidity |T|. The BCS wave function does have ODLRO 
via cooper pairing 0, however, it is not an eigenstate of the Hamiltonian. In 
the search for the mechanism of high temperature superconductivity, the Hub- 
bard model is studied intensively. Yang constructed explicitly many eigenstates, 
which is metastable for attractive interaction, of the Hubbard Hamiltonian possess- 
ing ODLRO [^]. Furthermore, Yang and Zhang uncovered an SO (4) symmetry and 
put forward a sufficient condition for any state to possess ODLRO [Q. Then Zhang 
predicted a triplet of collective modes, the massless Goldstone mode and a pair of 
massive modes, in the superconducting state ||T3f . 

In this paper we extend these discussions to the Hubbard bilayer, i.e., two Hub- 
bard planes coupled by interlayer hopping. The motivation is that strong coupling 
between the adjacent C11O2 planes of one double plane has been directly revealed 
experimently || 1§ . The rj pairing operator generalized from Yang's is intro- 
duced in Sec. II. Explicit eigenstates with ODLRO are constructed in Sec. III. 
SO (4) symmetry and the sufficient condition for a state to possess ODLRO is dis- 
cussed in Sec. IV, where a constraint relation among different kinds of ODLRO is 
also given. In Sec. V, the triplet of collective modes are discussed. A brief summary 
is contained in Sec. VI, where some open problems are raised. 
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II. rj pairing 

The Hamiltonian is 

H = T\\(A) + T\\{B) + T ± + V, (1) 

r,S 

= Y, e ( k )( a U a k] + a{ l a kl ), (2) 

k 

T,|(B) = t||£(&I T VH T + tfiW 

r,S 

= Y,< k )( b lM + b lM), (3) 
k 

T± = tj_ XX a rT fo rT + <Al b rl + ^T<M + ^j<M) 
r 

= t ± £(e-^4 T 6 fcT + e-^alM + e*"%a fcT + e*'%a kl ), (4) 

k 

V = U Y\( a U a r] ~ V)( a li a ri -fi) + {b^b r1 - fi){b j rl b rl - (j,)], (5) 

r 

where r is a two-dimensional integral coordinate variable on each L x L lattice, 5 is its 
nearest-neighbors on the plane, k is a two-dimensional integral momentum variable. 
q z denotes the momentum perpendicular to the planes, the distance between which 
is d = za — zb, where Za and zb are the vertical coordinates of the planes A and B, 
respectively. a r | and a r j are the coordinate-space annihilation operators for spin-up 
and spin-down electrons on layer A, respectively, b r ^ and b r ^ are those on layer B. 
a^, CLhii frfct and b k i are the corresponding annihilation operators in the momentum 
space. pJJ is the chemical potential. Tji (A) and Tji (B) are kinetic energies on the two 
layers, respectively. Tj_ is the interlayer coupling. V is the on-site electron-electron 
interaction potential energy. Following Ref. ||, we use 

e(k) = 4 — 2 cos k x — 2 cos k y (6) 



to make the Tji(A) and T\\(B) positive. Actually as has been pointed out ||, only 
e(k) + e(7r — k) — constant is required to make r\ pairing possible. 
For our present model, the operator 77 is introduced as 

rj = 7711(A) +m (B)+rj ± , (7) 

with 

r 

= e 2i ^ ZA Y. a ^-H, (8) 
k 



771, (B) = ^2^'^ 



rl 

= e**"° ^vK-ki, (9) 

k 

and 

r 

= e^ ZA+z ^J2( a klbn- kl -a kl b n ^), (10) 

k 

where 7r is two-dimensional. Actually 7711(A) and 7711 (B) are on-site 7r-momentum 
pairings, while 77^ is the interlayer nearest-neighbor 7r-momentum pairing. The 
definition is, of course, only meaningful when L is even. 
One can obtain the commutators 

^(A)+mB), m (A)+ m (B)) = -St,, [7711 (A) +77||(5)], (11) 



[Tn(A)+I||(S),77 ± ] = -8t m , (12) 
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[71, 7j||(A) + 7j,| (B)] = -2t ±V± , 

[71, 77J = -2t ± [ m (A)+ m (B)}. 
Under the constraint condition 

a^a r ^ + blfi r i = a>l^a r i + b\^b r ^ = 1, 

we may obtain 

[V, 7711(A) + 7711(B)] = -(l-2/i)C/(77|,(A)+ 7711(B)), 

[V,v±] = -(l-2fi)U V± . 
Therefore under the constraint condition flT5|), we have 

[h, v ] = -Er), 

where E = 8t\\ + 2t± + (1 - 2fi)U. 
The total momentum operator P is 

P = H( k + qz)(4t a fcT + a li a H + 4t & *T + b H b H), 

k 

the commutator of which with 77 is 

[P, 77] = -TFT/. 

III. Eigenstates with ODLRO 

We can generate the state 

i>N = f3(rf) N \vac >, 
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where \vac > is the vacuum state 



(3 = (N\(M - N)\/M\) 1/2 , M = L 2 (22) 



is the normalization factor. Because of (|I8D and (p0|), ipN is a simultaneous eigenstate 
of if and P, 

#ViV = NEipjy, (23) 
PViv = NirvpN- (24) 

Just similar to Ref. 0, one can construct another state ip' N , by replacing rj in 
Eq. ( |2T| ) with the 0-momentum pairing operator 

r?o = XX^fM + b r ]Ki + a r -\b r i — a r ib r i). (25) 

r 

V>at nas the same expectation value for H &s ipN but is not an eigenstate. So t/V 
is not the ground state. Similar to Ref. ||, it can be proved to be metastable if 
U < 0. 

Now we examine ODLRO. By definition, a state ip possess ODLRO for a local 
operator Q r means that the ODLRO correlation function 

lim i) ] Q ] s Qr^ ^ 0. (26) 

\r— s[— >oo 

In our system, there are three kinds of ODLRO correlation functions, the in- 
tralayer one for on-site pairing ^j v -atoLa r -[-o r j.'0jv = V'iv^si^st^T^r|'0Ar, the inter- 
layer one for on-site pairing ^jv^li &st a rT 'MV'iV = ^jv^^T^T^*!^' anc ^ that f° r 
interlayer pairing ^(fe^aL — 6Lat)(a r i-6 r j — a r jfe r |)^Ar, r ^ s. They all equal to 

M(M-l) ' 1 J 

Therefore three kinds of ODLRO are simultaneously possessed by ip^- 



IV. SO (4) symmetry 



Define \ = J + = J x + U y , \ = J_ = J x - iJ y , J z = \'Y^{a\ ] a k ^ + a j kl a kl + 

k 

4t^T + b kib kl + a! T 6fcT + a li b ki + 4t a fc*r + b ll a H) - 0ne ma y find tnat ^x, ^ 
and J z are generators of an SU(2) symmetry. Under the constraint condition Eq. 

(El), 

[H, J 2 } = [H, J z \ = 0. (28) 
If the total momentum operator is trivially re-defined as 

P' = $^(k - + q z )(4 T a fcT + a ] ki a kl + b\^b k} + 6^6 fci ), (29) 
fc 

we have 

[P', J] = 0. (30) 

From flUf ) we know that different from simple Hubbard model where [H, J] = 
[Q], we cannot have [H, J] = for coupled bilayer even if we re-define e{k) as 
— 2cos/c x . — 2 cos k y and let fi = 1/2, since [T±, rj\ cannot vanish. 
The particle-hole pairing operator is 

C = H(«rT«ri + Mr|) = HKt4| + M^)" ( 31 ) 
r fc 

Defining £t = J' x + iJ' y , ( = J' x - J y , J' z = ^Y,( a H a U + A^-a^a^ - &£ T & fcT ), 

k 

one may find that J' x , J' y , J' z are generators of an SU(2) symmetry, and that 

[H, J'] = [P', J'] = 0. (32) 

Note that the constraint condition Eq. ( pT5[) is not required for Eq. (0), which just 
represents the SU(2) symmetry of spin. 
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Consequently, as in the simple Hubbard, there is also an SO (4) symmetry in the 
Hubbard bilayer, with the generators properly defined. Most of the discussions in 
Ref. 0] can thus generalized in a straightforward way, for example, many eigenfunc- 
tions for H and P' can be obtained explicitely, the symmetry properties of energy 
spectrum are almost also valid. One only needs to take into account that now there 
are two planes and the system is also symmetric under the exchange of the two 
planes. 

However, the sufficient condition for a state to possess ODLRO should be re- 
considered. It can be stated as follows. 

Theorem. For any state ifi, if j 2 — j 2 = 0(M 2 ), where j and j z are the quantum 
numbers of J 2 and J z , respectively, then there is at least one kind of ODLRO 
among intralayer one for on-site pairing, interlayer one for on-site pairing and that 
for interlayer nearest-neighbor pairing. 

Proof. Assume matrix Q has matrix element which is the ODLRO correlation 
function of the operator A, 

Q sr = < A r \Q\A s > = < VIA^A^ >, (33) 

Using 

< A s \<f) >= M~ ll2 e™- S (34) 
as the trial wavefunction for Q, we obtain the expectation value of Q 

<Q> = <0|Q|0>=^£ eMs ~ r) < ^AjA^ > • (35) 
On the other hand, we have 

4(j 2 - f z + J + jz) = < V|f?^|V > = E e^ s ~ r) (2Q^ + 2Qi 2 r ) + QfJ), (36) 
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where 



Q { sr = < ^Wll a l-\ a ^arl\^ > = < Iplbl^b^br^lp > (37) 

Qf) = < i)\a} si a\^b rl \i) > = < i)\b ] si b\^a rl \i) > (38) 
Qf) = < 4)\{b ] sl a^ - b\^a si ){a r] b rl - a H b r ^ > . (39) 



Therefore 



2 < g« > +2 < > + < g( 3 ) > = l(j 2 - f z ) + 0(1). (40) 

If j 2 — il = 0(M 2 ), at least the largest of the eigenvalues of and is 

O(M), according to Yang's Theorem there is corresponding ODLRO. Q.E.D. 

Tian proved a theorem that the ODLRO correlation function of a local operator 
R r decays if there is another local operator S r satisfying [H, S r ] = R r f|. Using 
this theorem, we can obtain a useful result giving constraint on the ODLRO. By 
calculating [H, X)°rt°i]> we know that t||Q (nn) + t±Q (3) - U(2fi - 1)Q (1) decays, 

r 

here and are as defined above, Q( nn > is the ODLRO correlation function for 
intralayer nearest-neighbor s-wave pairing ^(a r -|-a r+( 5j — a r ^a r+ s^). So the existence 

r,S 

of one of the three ODLRO implies at least one of the other two. Furthermore, 
if fi = 1/2, interlayer pairing ODLRO exists if and only if ODLRO for intralayer 
nearest-neighbor s-wave pairing exists. Consistent result has been obtained for the 
pairing amplitudes in the states with spontaneous gauge symmetry breaking [ITU . 
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V. Collective modes 



In the Hubbard model, the SU(2) pseudospin symmetry contains the U(l) phase 
symmetry as a subgroup, whose spontaneously breaking, i.e., the superconductivity, 
gives rise to a pair of massive collective modes which together with the usual Gold- 
stone mode form a triplet representation of the pseudospin symmetry [T3j . Here we 
may see that these arguments can also extended to the bilayer model. 

First it can be found that the operators A_ = ^(a r |a r j + b r ^b r i + a r ^b r i — 

r 

a rl b r1 ), A+ = -AL, and A = \*Y^(a\ CJ a ra + b\ a b Ta + a\ a b ra + b\ a a ra ) form an 

rcr 

irreducible tensor of rank 1 under the SU(2) defined by J x , J y and J z in the last 
section. Then consider the response function 

D a (t,t') = -±9(t-t') < [J a (t),A a (t')} >, (41) 

where a denotes +, - or 0. One may obtain 

Do = (42) 
uj + id 

D±(uj) = (43) 
uj ± hj + io 

Therefore if the U(l) symmetry generated by J z is spontaneously broken, i.e., 
if there is superconductivity with intralayer or(and) interlayer pairings, < A + > = 
— < A_ >* 7^ 0, then there is a triplet of collective modes with energies ±E and 0, 
respectively. 
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VI. Discussions 



To summarize, by introducing an rj operator generalized from Yang's and under 
a constraint condition, the discussions concerning ODLRO and SO (4) symmetry 
for simple Hubbard model are extended to Hubbard bilayer. In the explicitly con- 
structed eigenstates, there are simultaneously three kinds of ODLRO, i.e., the in- 
tralayer one for on-site pairing, the interlayer one for on-site pairing, and that for 
interlayer nearest-neighbor pairing. By properly defining the generators, it is found 
that there is also an SO (4) symmetry. The sufficient condition for a state to possess 
ODLRO is that at least one of the three kinds of ODLRO exists if the condition 
is satisfied. We also obtain a constraint relation among the ODLRO for interlayer 
nearest-neighbor pairing, the intralayer ODLRO for on-site pairing, and that for 
intralayer s-wave nearest-neighbor pairing. 

Though the explicit eigenstates possess the three kinds of ODLRO simultane- 
ously, the sufficient condition can only ensure that at least one of them exsits. 
Another kind of ODLRO involves in the constraint relation. Further studies are 
needed to clarify whether these kinds of ODLRO must exist simultaneously. 

Some exact results about the ground state of the simple Hubbard model are 
based on the exchange of 77 (pseudospin) and £ (spin) operators under particle-hole 
transformation of the spin-up electrons [O] [I2~|. This situation is lost for the bilayer 
since the interlayer pairing is involved in our 77 operator. 

Recently the collective modes predicted by Zhang |TB|] was generalized to the 
triplet case, which is then claimed to be responsible for the 41 meV peak in the 
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neutron scattering spectrum on Y Ba<iCuO x f[4" |. However the dependence on q z 
U could not be explained. Our results suggest that the bilayer model might be 
incorporated into the framework of this theory to account for the q z dependence. 
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